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c/3 . Abstract 

New, gauge- independent, second order Lagrangian for the motion of classical, charged test 
particles is proposed. It differs from the standard, gauge-dependent, first order Lagrangian 
by boundary terms only. A new method of deriving equations of motion from field equations 
Qh| is developed. When applied to classical electrodynamics, this method enables us to ob- 

tain unambiguously the above, second order Lagrangian from the general energy-momentum 
conservation principle. 
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The motion of classical, charged test particles, in the classical Maxwell held is derived usually 



Q^ , from the gauge- dependent Lagrangian function 



CZJ - 

O | L = Lpartide + L int = -VI - v 2 (m - eu^A^t, q)) , (1) 

J>-^' where denotes the (normalized) four-velocity vector 

(un = (u°,u k ):=^==(l,v k ), (2) 

and v k := q k (we use the Heaviside-Lorentz system of units with the velocity of light c = 1). 

Since the Lorentz force eu^f^ u derived from this Lagrangian is perfectly gauge invariant, it 
is not clear, why we have to use the gauge- dependent interaction term eu^A^, with no direct 
physical interpretation. Moreover, in this approach equations of motion are not uniquely implied 
by held equations. As an example, non-linear forces of the type u^f^u^xf^ cannot be a priori 
excluded. 
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In the present paper we show that the same Lorentz force may be derived from a gauge- 
invariant, second order Lagrangian C: 



£ = L partide + Cua = -Vl-v 2 (m - aWMfiit, q, v)) , (3) 

where a 11 := u u V u u^ is the particle's acceleration. The skew-symmetric tensor M*™ (i, q, v) is 
equal to the amount of the angular-momentum of the field, which is acquired by our physical sys- 
tem, when the (boosted) Coulomb field ff%, , accompanying the particle moving with constant 
velocity it through the space-time point y = (t, q), is added to the background (external) field. 
More precisely: the total energy-momentum tensor corresponding to the sum of the background 
field and the above Coulomb field decomposes in a natural way into a sum of 1) terms 
quadratic in f^ u , 2) terms quadratic in tj^,' u ^ and 3) mixed terms. The quantity MJ%f is equal to 
this part of the total angular- momentum M^ u , which we obtain integrating only the mixed terms 
of the energy-momentum tensor. It is proportional to the particle's electric charge e contained 
in the particle's Coulomb field. The interaction Lag rangian Hint can thus be expressed in terms 
of the following quantity 

Qti ■= \u"M$ , (4) 

which is obviously orthogonal to the four- velocity u^. Therefore, it has only three independent 
components. In the particle's rest frame we have Q$ = and the three-vector eQk equals to the 
amount of the static moment acquired by the system when the particle's own field is added to 
the background (external) field. Hence, we have from (0) and (0): 



C m t = Vl-v 2 ea^Q^t, q, v) . (5) 

We prove in Section [2] that, for a given external field /^„, the new interaction Lagrangian 
(H) differs from the old one by (gauge-dependent) boundary corrections only. Hence, both 
Lagrangians generate the same physical theory (although the new Lagrangian is of second dif- 
ferential order, its dependence upon the second derivatives is linear; this implies that the corre- 
sponding Euler - Lagrange equations of motion are of the second order). The relation between L 
and C is, therefore, analogous to the one well known in General Relativity: the gauge- invariant, 
second order Hilbert Lagrangian for Einstein equations may be obtained starting from the first 
order, gauge-dependent Lagrangian and supplementing it by an appropriate boundary term. 

At this point, our result can be summarized as follows: a physical interpretation of the inter- 
action Lagrangian eu^A^ has been found. Up to boundary terms, it is equal to the interaction- 
angular-momentum a^u^M™ 1 . The question arises whether such an identity is a pure coincidence 
or is a result of a universal law of physics. 

In the second part of the paper (Sections 3 - 5) we try to convince the reader that the second 
conjecture is correct. In fact, we propose a new method of deriving equations of motion from 
field equations. The method is based on an analysis of the geometric structure of generators of 
the Poincare group, related with any special-relativistic, lagrangian field theory. This analysis 
leads us to a simple theorem, which we call "variational principle for an observer" (see Section 
^). Applying this observation to the specific case of classical electrodynamics, we show how 
to derive, in principle, equations of motion from field equations. This derivation is based on a 
following heuristic idea: a point-like particle has to be understood as an approximation of an 
extended particle, i. e. of an exact, stable, soliton-like solution of a hypothetical theory of matter 
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fields interacting with electromagnetic field. To prove that the Maxwell theory (describing the 
free field outside of the "soliton's strong-field-core" ) implies equations of motion for the solitons 
themselves, we need several qualitative assumptions of heuristic nature, about their stability. 
Under these assumptions the gauge-invariant Lagrangian describing the motion of point particles 
is unambiguously derived from the general invariance principles of the theory. 

The mathematical status of the above derivation is, therefore, similar to the Einstein - Infeld 
- Hoffmann derivation of equations of motion from field equations in General Relativity. It does 
not depend upon a particular model which we take for the description of the gravitating body 
under consideration (e. g. a hydro dynamical or an elastomechanical model). The derivation 
is valid for any stable body and enables us to describe (in a good approximation) its motion 
in a model-independent way, as a geodesic motion of a point-particle. Hence, even if we have 
at the moment no realistic mathematical theory describing the interior of a star and fulfilling 
all the necessary assumptions, we can expect that the above equations of motion are valid for 
relatively stable objects. The present paper shows that a similar argumentation is possible also 
in electrodynamics: Lorentz force acting on test particles does not need to be postulated as an 
independent physical law, but is implied by the geometry of Maxwell field, provided one accepts 
the existence of the hypothetical fundamental theory of matter fields, admitting sufficiently 
stable soliton-like solutions. 

The above result follows immediately from the consistent theory of interacting particles and 
fields (cf. p|], P|), called Electrodynamics of Moving Particles. All the formulae of the present 
paper can be derived directly from the above theory in the test particle limit (i.e. m — > 0, e — > 
with the ratio e/m being fixed). The present paper, however, does not rely on this, much more 
general, context. The consistent theory of test particles is constructed independently. 

The paper is organized as follows. Section |2| contains a direct proof that our new Lagrangian 
differs from the standard, gauge-dependent one by boundary terms only. Section |3| contains our 
basic geometric observation concerning any relativistic, Lagrangian field theory, which makes 
our approach possible. Using it, we give in Sections || and || the derivation of our variational 
principle. In particular, the renormalization procedure defined in Section || depends upon the 
stability assumptions (which may be relatively difficult to check for any specific mathematical 
model) . Accepting these assumptions on a heuristic level we obtain our theory of point particles 
uniquely, as a necessary consequence of energy-momentum conservation. 

In a forthcoming paper we are going to present the gauge-invariant Hamiltonian structure 
of the above theory. 

2 Equivalence between the two variational principles 

The easiest way to prove the equivalence consists in rewriting the field fav in terms of the 
electric and the magnetic induction, using a special accelerated reference system, adapted to 
the particle's trajectory £, which we define in the sequel. The system will be also very useful to 
formulate our "variational principle for an observer" in the next Section. 

We begin with any laboratory coordinate system (y^) = (y°,y fc ) in Minkowski space-time 
M and parameterize the trajectory £ with the laboratory time t = y°. Let y k = q(t) be 
the corresponding coordinate description of £. At each point of the trajectory we choose an 
orthonormal tetrad (e^), such that its element eo is tangent to £, i. e. is equal to the four-velocity 
vector u. Take now the unique boost transformation relating the laboratory time axis d/dy° 
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with the observer's time axis eo- We define the vector by transforming the corresponding 
d/dy k - axis of the laboratory frame by the same boost. 

It is easy to check (cf. [|J) that the above definition implies the following, explicit formula: 

_i_| 5 + (4+ , (vV „ t) | r , (6) 



where we denote tp(z) := \ ( ^\_ z — l) 



i 



/T^I(l+v / T 3 i) ' 

Finally, we parameterize space-time points by four coordinates (t,x k ): 

y(t,x):=(t,c l (t))+x k e k (t) . (7) 

Using (||) we obtain the following relation between our curvilinear parameters (t,x k ) and the 
laboratory (Lorentzian) coordinates (y^): 

y°(t,x l ) := t+ 1 x l vi(t) , 
V 1 - v W 

y k (t,x l ) := q k (t) + (d k + i p(v 2 )v k v l )x l . (8) 

The above formula may be used as a starting point of the entire proof. To rewrite field equa- 
tions with respect to this system we need to calculate the components of the flat Minkowskian 
metric in our new coordinates. We see from (^) that, for a given value t, parameters (x k ) are 
cartesian coordinates on the 3-dimensional hyperplane Sj, orthogonal to £ at the point (t, q(i)). 
Hence, we get g k i = $ki for the space-space components of the metric. 

The remaining information about the metric is carried by the lapse function and the shift 
vector, which again may be easily calculated from formula (||): 

N = _L_ = Vl-v2 (l + OiX*), (9) 

V-g m 

N m = gom = Vl - v 2 e mH uj k x l . (10) 

Here, by a 1 we denote the rest-frame components of the particle's acceleration. They are given 
by formula: 

— e = — u = a l ei , (11) 
ar ar 



si + ^yvk)^ (12) 



where r is the proper time along (. A straighforward calculation gives us (formulae (J2|) and (|60 
the following value: 

i 1 
a =1^ 

where v k is the acceleration in the laboratory frame. 

Moreover, at each point of £ we define the rotation vector (u 3 ) of the tetrad by the following 
formula: 

= meo - e^tJek . (13) 
Again, straighforward calculation leads to the following expression: 

uj m = 1 <p{~v 2 )v k v l e k i m . (14) 
V 1 — v z 
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The transformation (||) is not invertible. Coordinates (x k ) are regular parameters on each 
Hf But hyperplanes corresponding to different values of t may intersect. Hence, the same 
physical event may correspond to different values of coordinates (t,x k ). 

Nevertheless, we may describe the free Maxwell field in terms of our parameters (t,x k ). In 
particular, we have: 

A k - d k A =: f 0k = -ND k + e mH N m B l , (15) 

where D k and B k are the electric and the magnetic field on each Ej, N and N m are given by 
(H) and ( |io|) (for the description of the Maxwell field with respect to an accelerated reference 
system see e. g. Q). 

Let us multiply fllq ) by ^ and integrate this scalar product over S t with respect to the 
Lebesgue measure dx 3 . We first integrate over the exterior of the sphere S(ro). Observe, that 
in this region the following identity holds: 

^d k A = d k ^ . (16) 
Moreover, we have: 

t,e mkl N m B l = -VT^ d l ( B 1 ^] . (17) 



r 

Hence, after integration, we obtain in the limit ro — ► 0: 

k k 

^A k d 3 x + A7rA (t,0) = - / ^ND k d 3 x (18) 

(Ao(t,0) is the only surface term which survives in the limit, due to the standard asymptotic 
behaviour of the field). Observe that the constant part of the lapse function (|9|) does not 
produce any contribution to the right-hand side of the above formula, because the flux of the 
field D k through any sphere S(r) vanishes due to the Gauss law. Hence, we may replace "iV" 
by "Vl — v 2 diX 1 " under the integral and obtain 

„k 



[ ^ND k d 3 x = —y/T^a k [ x k B n D n d 3 x 
JEi r 6 e 7s t 



4-7T 4-7T 

.—a k Mi^ = -—a^M™ 1 , (19) 



where 



D n := (20) 



e x 
4-7T r 

is the Coulomb field on St, corresponding to the charge e. 

The lower index "0" in our particular system comes from the vector which is proportional 
to the particle's velocity u^, where the proportionality coefficient Vl — v 2 is due to the ratio 
between the proper time and the laboratory time on the trajectory. This means that Ao(t, 0), 
calculated in our particular coordinate system, is equal to Vl — v 2 u^A^t, q(i)) in any other 
coordinate system. The same is true for M 1 Jq . We have, therefore: 



e 

47T 



/ t,A k d 3 x + VT^ 2 eu^AJt, q(i)) = VT 3 ^ aV 'M™ 1 . (21) 
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Integrating this identity over a time interval [^,£2] we finally obtain 

f £ - = £ Lmt - h (k i Ak ** " 1, ^ - 41 A ) ' (22) 

i. e. both Lagrangians differ by boundary terms only. 

We will see in the sequel that the value of M^™* does not depend upon the choice of a 
particular hypersurface Xf which we have used for integration. Any other £, which intersects 
the trajectory at the same point and is flat at infinity will give the same result. We conclude that 
there is always a gauge in which both C and L coincide, e. g. the gauge in which the monopole 

k 

part of the radial component A r := A^^- vanishes. 



3 Variational principle for an observer 

Consider any relativistic- invariant, Lagrangian field theory (in this paper we will consider mainly 
Maxwell electrodynamics, but the construction given in the present Section may be applied to 
any scalar, spinor, tensor or even more general field theory). Choose any non-inertial observer, 
moving along a time-like trajectory We want to describe the field evolution with respect to 
the observer's rest frame. For this purpose we choose the space-time parameterization defined 
in the previous Section. 

The field evolution with respect to the above non-inertial reference frame is a superposition 
of the following three transformations: 

• time-translation in the direction of the local time-axis of the observer, 

• boost in the direction of the acceleration a k of the observer, 

• purely spatial 0(3)-rotation tu m . 

It is, therefore, obvious that the field-theoretical generator of this evolution is equal to 

H = y/1- v 2 (£ + a k R k - uj m S m ) , (23) 

where £ is the rest-frame field energy, R^ is the rest-frame static moment and S m is the rest- 
frame angular momentum. The factor y/l — v 2 in front of the generator is necessary, because 
the time t = x°, which we used to parameterize the observer's trajectory, is not the proper 
time along £ but the laboratory time. For any point (t, q(t)) € C the values of the Poincare 
generators £,Rk and S m are given as integrals of appropriate components of the field energy- 
momentum tensor over any space-like Cauchy surface £ which intersects £ precisely at (i,q(i)) 
(due to Noether's theorem, the integrals are independent upon the choice of such a surface). 
These values are, therefore, equal to the components of the total four-momentum p^ and the 
total angular momentum M^ u of the field, calculated in the observer's co-moving frame e M , 
i.e. £ = —p^, R^ = —M fJiU u u and = 0. 

Given a field configuration, we are going to use the quantity if as a second order Lagrangian 
for the observer's trajectory. For this purpose let us first choose a "reference trajectory" £o- 
Next, for each point (i,q(i)) G Co let us calculate the corresponding "reference values" of the 
generators £(t), Rk(t), S m (t) and insert them into H. Finally, consider the function ([23| ) obtained 
this way as a Lagrangian depending upon a generic trajectory £ via its velocity v and acceleration 



v, according to (12) and (14). 
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Theorem 

Euler- Lagrange equations derived from the above Lagrangian are automatically satisfied by the 
trajectory C = Co- 

This theorem was derived in M in a much more general framework. Within this framework, 
it was an obvious consequence of the invariance of the theory with respect to the choice of an 
observer. More precisely, the function "—if" was proved to be a Routhian function playing the 
role of the Lagrangian with respect to the observer's degrees of freedom and the Hamiltonian 
(with opposite sign) with respect to the field degrees of freedom. For purposes of the present 
paper we do not need, however, this larger context. The Theorem may be checked by simple 
inspection: Euler-Lagrange equations derived from the second order Lagrangian (|23| ) are au- 
tomatically satisfied as a consequence of the field energy-momentum and angular-momentum 
conservation (this direct proof was also given in |@|). 



4 Adding a test particle to the field 

From now on we limit ourselves to the case of electrodynamics. This means that the field energy- 
momentum and angular-momentum are defined as appropriate integrals of the components of 
the Maxwell energy-momentum tensor 

n = r x u - \wr x f KX . (24) 

Suppose now that to a given background field f^ we add a test particle carrying an electric 

charge e. Denote by fj$J the (boosted) Coulomb field accompanying the particle moving with 
constant four-velocity u, which passes through the space-time point y. Being bi-linear in fields, 
the energy-momentum tensor T total of the total field 

C al -'=f^ + ^ U) ( 25 ) 

may be decomposed into three terms: the energy-momentum tensor of the background field 
r /ieid the Coulomb 

energy- momentum tensor x partlcle , which is composed of terms quadratic 
in and the "interaction tensor" T mt , containing mixed terms: 

rptotal rpfield _|_ rpparticle _|_ rpint (26) 

Unfortunately, T total does not lead any longer to a globally conserved quantity, as in the case 
of a relativistic invariant field theory. Indeed, the relativistic invariance has been broken by the 
choice of the electric current localized on the trajectory. Moreover, due to the Coulomb field's 
singularity, integrals which would have been necessary to obtain the generator (23) from (|26|) 
are ill defined. Hence, techniques introduced in Section 3 cannot be used directly. 



For reasons which will be fully explained in the next Section, we replace (23) by the following, 
well defined, "renormalized" quantity. To obtain this quantity we will integrate first two terms 
of ( |26| ) over any £ which passes through y = (t, q(t)). Integration of T^ ield (no singularity) and 
T mt (an r~ 2 - singularity) is possible and gives well defined quantities, which we call respectively 
II field ^ "background held generator" ) and H int ( "interaction generator" ) . Because the left-hand 



side and the first two terms of the right-hand side of (26) have a vanishing divergence (outside 
of the particle's trajectory), we conclude that also T has a vanishing divergence. This implies, 
that the above integrals are invariant with respect to changes of E, provided the intersection 
point with the trajectory does not change (see |Q] for more details). 
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Unfortunately, the Coulomb tensor T parttcle has an r -4 singularity at y and cannot be in- 
tegrated. According to the renormalization procedure defined in [|j and sketched briefly in 
the next Section, we replace its integrals by the corresponding components of the total four- 
momentum of the particle: pP^ rUcle = mu\ and the total angular momentum: ]\/[P% rtlde = 0. We 
define, therefore, the renormalized particle generator as follows: 

HPZ Ude = Wl-v 2 . (27) 
Consequently, we define the total (already renormalized) generator as a sum of three terms 

jjtotal jjfield _|_ jjint _|_ jjparticle (28) 

where the first term is quadratic and the second term is linear with respect to the background 
field f^. 

Let us observe that the only non-vanishing term in H mt comes from the static moment 
term R in (|23|), because the mixed terms in both the energy £ and the angular momentum S 
vanish when intergated over any E. The easiest way to prove this fact consists in choosing the 
hypersurface £t which is orthogonal to the velocity u at (i, q(i)), i.e. the rest-frame surface 
(our integrals do not depend upon the choice of a hypersurface). On this surface, the Coulomb 
field fjn; ^ is spherically symmetric and carries, therefore, only the monopole component. On 
the other hand, the monopole component of the background field fun vanishes as a consequence 
of the homogeneous Maxwell equations (no charges!). The mixed term in the energy integral 
is, therefore, a product of a monopole-free functions and the pure monopole. Hence, it vanishes 
after integration. A similar argument applies to the angular momentum S. 

Finally, we have defined 

h%£> = H* ieid - a ^Mjz\t, q , v) + yr^2 m ( (29) 

where the interaction term is defined as the following integral 

Mfiiv) ■■= I y,)Ttf(x) - (x v - yu)T^(x)} dS A (x), (30) 

J S 

and S is any hypersurface which intersects the trajectory at the point y = (t, q(t)). 

In particular, using the particle's rest-frame and integrating over the rest-frame hypersurface 
£f we obtain formula (|l9|) for Mjff. 

As we have already mentioned, the quantity i?*°n defined this way, cannot be directly used 
in the framework of pure electrodynamics, as the "observer's Lagrangian". In the next Section 
we will show, however, that this quantity provides a good approximation of the total generator 
H within a more general framework. In this framework particles are no longer point-like, but 
are extended objects, described by matter fields interacting with electrodynamics. We will 
conclude, that within this approximation our "renormalized generator" may be used to derive 
an "approximative trajectory" of such extended objects. 



5 Renormalization. Derivation of equations of motion from field 
equations 

According to the approach developed in M and [g] we treat the moving particle as a solution of 
a hypothetical "fundamental theory of matter fields interacting with electromagnetic field" . We 
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assume that such a theory is a relativistic, Lagrangian (possibly highly non-linear) field theory. 
Moreover, we assume linear Maxwell theory as a limiting case of the above theory, corresponding 
to sufficiently weak electromagnetic fields and vanishing matter fields. 

We will suppose that the particle, whose interaction with the electromagnetic field we are 
going to analyze, is a global solution of the above field theory, having following qualitative 
properties: 

1. it contains a tiny "strong field region" , concentrated in the vicinity of a time-like trajectory 
C, which we may call an approximate trajectory of the extended particle, 

2. outside of this strong field region the matter fields vanish (or almost vanish in the sense, 
that the following approximation remains valid) and the electromagnetic field is sufficiently 
weak to be described by linear Maxwell equations. 

To be more precise, we imagine the "particle at rest" as a stable, static, soliton-like solution 
of our hypothetical "super theory". The solution is characterized by two parameters: its total 
charge e and its total energy m. The energy is not concentrated within the interior of the particle 
but contains also the part of the energy carried by its "Coulomb tail" . This means that m is 
an already renormalized mass, (or dressed mass), including the energy of the field surrounding 
the particle. Within this framework questions like "how big the bare mass of the particle is and 
which part of the mass is provided by the purely electromagnetic energy?" are meaningless. In 
the strong field region (i. e. inside the particle) the energy density may be highly non-linear and 
there is probably no way to divide it consistently into two such components. 

Due to relativistic invariance of the theory, there is a 6 parameter family of the "uniformly 
moving particle" solutions obtained from our soliton via Poincare transformations. 

Now, an arbitrarily moving particle is understood as a "perturbed soliton" . This means that 
it is again an exact solution of the same "super theory" , with its strong-field-region concentrated 
in the vicinity of a time-like world line which is no longer a straight line, as it was for "uniformly 
moving particles". Let us choose an observer who follows this "approximate trajectory" £. We 
know that he automatically satisfies the Euler-Lagrange equations derived from the second order 
Lagrangian (p^), where £ , R k and S m are the quantities calculated for the complete non- linear 



Suppose now, that the particle may be treated as a test particle. This means, that the total 
field outside of the particle does not differ considerably from a background field /, satisfying 
homogeneous Maxwell equations. Using this hypothesis we may approximate, for each point 



Here, by we denote the total energy-momentum of the "super theory", corresponding to 
the "uniformly moving particle" solution, which matches on £ the position and the velocity 
of our particle. Stability of the soliton means that the the "moving particle solution" does not 
differ considerably from the "uniformly moving particle solution" inside of the strong field region 
(i. e. inside the particles). Hence, the contribution of the first term (T — T) to the integrals £, 
R k and S m may be neglected "inside the particle", i. e. we may replace it under integration by 
the purely Maxwellian quantity 



theory. 




(V 




rint 



(32) 
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As a result, we obtain the first two terms of (|29|). On the other hand, integrating the last term 
T in (31) we obtain without any approximation the corresponding value of the four-momentum 
of the particle. This way we reproduce the last term of (pS|), 

Replacing ( p3|) by its approximate value (^9|) and using the Theorem we conclude that the 
trajectory £ has to fulfill Euler-Lagrange equations derived from (|29|). Applying the Theorem to 
the linear Maxwell theory we conclude, that the term Hf ield produces Euler-Lagrange equations 
which are automatically fulfilled by (. Hence, we may drop out this term, leaving only the 
remaining two terms. They finally give us our formula (||) for the Lagrangian of the test particle 
(the sign has to be changed because as a Lagrangian we should have taken U —H" instead of 
U H" - see remark at the end of Section 0). 
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